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ABSTRACT 


Radar signal processing applications frequently require an estimate of 
the Doppler centroid of a received signal. The Doppler centroid estimate is 
required for synthetic-aperture radar (SAR) processing. It is also required 
for some applications involving target-motion estimation and antenna pointing 
direction estimation. In some cases f the Doppler centroid can be accurately 
estimated based on available information regarding the terrain topography, the 
relative motion between the sensor and the terrain, and the antenna pointing 
direction. Often, the accuracy of the Doppler centroid estimate can be 
improved by analyzing the characteristics of the received SAR signal. This 
kind of signal processing is also referred to as clutterlock processing. This 
publication reports on a Doppler centroid estimation (DCE) algorithm which 
contains a linear estimator optimized for the type of terrain surface that can 
be modeled by a quasi-homogeneous source (QHS) . 

Information on the following topics is presented in this publication: 

• an introduction to the theory of Doppler centroid estimation 

• analysis of the performance characteristics of previously reported 
DCE algorithms 

• comparison of these analysis results with experimental results 

• a description and performance analysis of a Doppler centroid 
estimator which is optimized for a QHS 

• comparison of the performance of the optimal QHS Doppler centroid 
estimator with that of previously reported methods 
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SECTION I 


INTRODUCTION 

The quality of SAR images strongly depends on the accuracy of the 
knowledge of the range history of the sensor relative to the terrain during 
the interval of observation. In many applications, this range (or signal 
phase) history can be accurately modeled by a second-order Taylor series 
expanded about the center of the synthetic aperture. The coefficient of the 
quadratic term is equal to the Doppler frequency rate. An error in this 
parameter causes broadening of the SAR impulse response. The coefficient of 
the linear term is equal to the Doppler frequency associated with a point 
located at the center of the antenna beam. An error in this parameter causes 
degradation of several image-quality factors including the signal-to-noise 
ratio (SNR) and the signal-to-ambiguity ratio (STAR) [1] . 

A SAR antenna pattern is usually symmetric with its peak gain at the 
center of the beam. The best SNR and STAR are obtained by selecting a 
processing-frequency band centered at the Doppler centroid. These image- 
performance factors are degraded when there is a mismatch between the center 
of the processing band and the Doppler centroid. In many other related 
applications, such as determination of the speed of moving targets [2, 3] and 
determination of antenna pointing direction [41, the system performance is 
directly affected by the accuracy of the Doppler centroid estimate. The value 
of developing a more accurate DCE algorithm is therefore apparent. 

Doppler parameters are often computed from the sensor-target relative 
position, velocity, and acceleration vectors [5] . Alternatively, these 
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parameters can be estimated from the SAR image data, as demonstrated by Li et 
al. [1]. The principle of this algorithm is that the intensity profile of the 
range-averaged power spectrum of the image data is a function of the antenna 
pattern in the azimuth dimension. The peak in the spectrum corresponds to the 
Doppler centroid. That method was found to enable an improvement in the 
Doppler centroid estimate when the ephemeris data contains significant errors. 
Such was the case in both the Seasat and SIR-B missions. Several other 
methods, similar in principle, were presented by Curlander et al. [6] . 

However, no performance analysis of those algorithms has been published. 
Consequently, comparison of performances among these algorithms in the 
available literature lacks a theoretical basis. 

In the first part of this publication, two previously reported DCE 
algorithms and a generalized algorithm are described. The performance 
analysis for this generalized algorithm is then presented. This analysis 
assumes that the statistical properties of the backscatter coefficients of 
targets follow those of a QHS . This assumption is valid for many physical 
sources. Hence, the result of this analysis is applicable to a large 
percentage of the SAR data. Experimental results obtained using Seasat and 
SIR-B data agree well with the analysis. 

One fact observed from the result of the performance analysis with a QHS 
is that the error of the Doppler centroid estimate provided by the generalized 
algorithm increases when the contrast within the SAR image increases. This 
indicates that the generalized algorithm is not optimized for the QHS. 

Further studies showed that this algorithm is also not an optimal approach for 
a strictly homogeneous source (HS) . This finding was the motivation for 
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developing an optimal DCE algorithm applicable to a wider range of terrain 
characteristics . 

The second part of this publication describes the design of an optimal 
DCE process. In this design, a number of optimal Doppler centroid measures 
are obtained by processing a number of spatially diversed subsets of image 
pixels with a linear estimator; each subset is strictly homogeneous. A single 
final Doppler centroid estimate is then given by the average of the weighted 
Doppler centroid measures. These weight coefficients are generated based on 
the maximum likelihood criterion so as to minimize the variance of the final 
Doppler centroid estimate. In the special case of the homogeneous source, 
these weight coefficients are proportional to the number of pixels in each 
subset . 

To demonstrate the validity of this optimal DCE algorithm and the 
improvement of the estimation accuracy, several Seasat and SIR-B SAR data sets 
were processed. Results show that the optimal DCE algorithm generally yields 
more accurate Doppler centroid estimates than the previous algorithms. 
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SECTION II 


DOPPLER CENTROID ESTIMATION FROM SAR DATA 

To illustrate the concept of the DCE by using SAR data, we will start 
with a simple case in which the SAR echo response is due to a single point- 
target (consisting of one scatterer with its dimension much smaller than the 
SAR resolution) . This concept is then extended to a more general case in 
which the SAR echo response is contributed by a large number of distributed 
targets (consisting of many scatterers in one resolution cell) . Previous DCE 
algorithms are briefly reviewed, and a generalized algorithm is then 
presented. 

A. POINT-TARGET SAR ECHO RESPONSE 

The echo response of a SAR from an isotropic point-target reflector can 
be modeled by its amplitude and phase variations in the along-track dimension 
[5], i.e., 

Rp (t ) = W y 2 (t) • exp{-j()>(t) } d) 

where W = (t) is the azimuth antenna weighting and $(t) is the phase history, 
which is a function of the sensor-target range history. The phase history can 
be modeled by a quadratic function, i.e., 

♦ <t> = 2rt ( f dc • t + | • f r • t 2 ) (2) 
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where f dc is the Doppler centroid and f r is the Doppler frequency rate. For 
this simple case, the Doppler centroid can be determined directly from the 
echoes by tracing the phase-change rate at the time of the response peak. For 
a response with a large time-bandwidth product , the amplitude of the Fourier 
transform of this response takes the same form. Therefore, one can also 
determine the Doppler centroid in the frequency domain. The echo response and 
spectrum in the azimuth dimension for a point-target are illustrated in 
Figure 1. Note that to satisfy Eq. (1), isotropy is not an essential property 
of the target. An adequate and less stringent condition is that the target 
backscatter be independent of the azimuth angle change over the illumination 
period. It is believed that many types of targets exhibit this property, 
especially for cases with a narrow antenna beamwidth. This property will be 
assumed in the following analysis . 


B. DISTRIBUTED-TARGET RESPONSE 


In determining the Doppler centroid for a distributed target, two 
problems may be encountered. First, tracing the phase change for each target 
by analyzing the response data in the time domain is impossible because the 
SAR echo response is a superposition of many point-target responses . Second, 
fading may exist in each spectral component because of the random-phase 
process of the naturally occurring targets. Consequently, the Doppler 
centroid estimate processed from a single azimuth power spectrum exhibits 
significant error. 
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Figure 1. SAR Response of a Point-Target: (a) Point-Target Echo Response 
in Azimuth, (b) Spectrum of Point-Target Echo Response in 




A commonly adopted method [7] is to obtain a few azimuth power spectra 


from the response data, average them to improve the signal-to-noise ratio, and 
extract the Doppler centroid by searching for the peak of the averaged 
spectrum. Li [1] noted that this method usually cannot provide high accuracy 
and showed that a much more accurate estimate can be obtained from the range- 
averaged spectrum of a processed SAR image. In the following, we summarize 
processing procedures of algorithms given by both Li et al. and Curlander et 
al . 


1. Algorithm I. (Li et al.) 

Process the SAR echo response into a full-resolution image with an 
initially estimated Doppler centroid and an accurate Doppler frequency rate. 
Fourier transform this image data in the azimuth dimension and perform 
detection to obtain the power density for each frequency channel. Average 
these detected spectra along the range dimension and search for a frequency 
channel which has the energy on both sides balanced. Iteration of this 
process may be required if the initial estimate used for processing the image 
contains a large error. 

2. Algorithm II. (Curlander et al.) 

Process the SAR echo response into four independent single-look images. 
Each image uses one-fourth the total processing bandwidth, and the spectra of 
these images are consecutive. Then compute the following quantity: 


2-4 



4 


where E i denotes the total energy in the i-th look. Repeat this procedure for 
several selected values of the Doppler centroid. Perform linear regression 
for E as a fuction of the processing Doppler values. The Doppler centroid 
estimate is selected from the intersection of the linearly fitted line and the 
line of AE = 0 . 

Both of these algorithms are illustrated in Figure 2. Algorithm I 
described above requires processing full-resolution SAR images, which is a 
disadvantage if the normal processor output is multi-look (reduced-resolution) 
imagery. Algorithm II involves a relatively large computation load and may 
not be efficient. Algorithm III (described below) is a generalized DCE 
algorithm which was evolved from the above algorithms. The performance 
analysis of algorithm III can be easily made and is applicable to the above 
algorithms . 

3. Algorithm III. 

The first step of this algorithm is to process the SAR echo response 
into N independent single-look images, each with 1/N of the total processing 
bandwidth, where N is an even number and the spectra of these images are 
consecutive. Then the following quantity is computed: 
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AVERAGED SPECTRUM OF A FEW 
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N/2 N 

X E i - I E i 

i=l i=N/2+l 

AE = (4) 

N 

i=l 

where denotes the total energy in the i-th look. Obtain the Doppler 

centroid estimate by subtracting the estimated error Af of the processing 
Doppler (Af - AE, is given in Eq. (35)) from the processing Doppler. 

Repeat this procedure until AE approaches zero. 

It is obvious that this algorithm is equivalent to algorithm I when N is 
equal to two. With N equal to four, this algorithm is very close to 
algorithm II. However, there is a slight difference in the performances, 
which will be discussed later. 
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SECTION III 


PERFORMANCE ANALYSIS OF DCE ALGORITHM III 

In DCE algorithm III, the accuracy of the estimated Doppler centroid is 
determined by the statistical characteristics of the normalized energy 
difference AE . It is therefore necessary to formulate the probability density 
function of AE . This pdf is used to derive the standard deviation of the 
Doppler centroid estimate. 

A. SAR IMPULSE RESPONSES OF INDEPENDENT LOOKS 


The SAR impulse response can be expressed as the convolution of the 
returned point-target echo with a reference function. If there is no 
amplitude weighting included in the reference function, the impulse response 
in the azimuth dimension is given by 

h(t) = Rp(t) @ {exp [ j4> (t) ] Rect(t/T)} (5) 

where (x) denotes the convolution, Rp(t) is the SAR point-target response 
after range correlation, and Rect ( . ) is the rectangular function defined by 

f 1 I t I < T/2 (fi v 

Rect (t/T) - ^ 0 otherwise 


where T denotes the processing aperture width. For a mismatch between the 
processing Doppler centroid and the true Doppler centroid by Af, the reference 
function is equivalent to the previous one having a time delay Af/f r , i.e.. 
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(7) 


h(t,Af) = Rp(t) (x) { exp [ j<|>(t) ] 

Rect ( (t - Af/f r ) /t) } 

For a reference function with a large time-bandwidth product (TBP) , the 
magnitude of the impulse spectrum takes the same form as the antenna weighting 
in the time domain, i.e., 

| F {h (t, Af ) } | - W a 1/2 ( (f - f dc ) /f r ) 

Rect ( (f - f dc - Af ) /PBW ) (8) 

where PBW is the processing bandwidth given by f r • T. It is convenient to 
define 


w a (f) = w a ( (f - f dc ) /f r ) 


(9) 


Images produced by a coherent system are corrupted by speckle noise. 
Methods of noncoherent averaging that reduce the speckle noise and are based 
upon properties of spatial or frequency diversity [8] have been reported 
[9,10]. The trade-off of speckle reduction is the broadening of the 
resolution width. In SAR, statistically independent looks can be generated 
from spatially separated subapertures. For the case with the total number of 
looks being equal to N, the impulse response corresponding to the i-th 
subaperture can be expressed by 


h i (t,Af) = Rp(t) @ { exp[ j<J) (t ) ] 

Rect ( (t - ti - Af / f r ) / (T/N) ) } (10) 
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where t i is the time delay associated with the i-th subaperture and 


t L = (i - 1) T/N - T/2 (11) 


If the TBP of each subaperture is also large, the spectrum of the i-th 
impulse response can be approximated by 

| F {h ± (t, Af ) } | = W a 1 / 2 (f) 

Rect ( (f - f L - Af ) / (PBW/N) ) (12) 


where 


fi = fdc + f r 


* ti 


(13) 


The following property can be derived directly from this approximation: 


hi (t, Af ) @ hj (t, Af ) = 0 for i ^ j 


(14) 


If the azimuth ambiguity (spectral aliasing due to limited pulse repetition 
frequency) makes a negligible contribution to the impulse response, the 
following properties can be shown to be true: 


E 


/• 

r 

| hi (t , Af ) 

2 dt = 

•j 


PBW 

2 


+ Af 


-PBW 

2 


+ Af 


w a (f) df 


(15) 
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and 


i 


r 


h £ (t , Af ) @ hj^* (t , Af ) 


dt 


/* + Af 

2 


W|(f) df 


— PBW 


+ Af 


(16) 


where * denotes the complex conjugate. 


B . TARGET MODEL 


The concept of the quasi-homogeneous source has been brought forth by 
Carter and Wolf [11] to model many physical light sources and to derive 
coherence properties and radiometric properties of light generated by such 
sources . 

In SAR, a similar model was used by Raney [12] to derive the SAR system- 
transfer function and some useful properties. With this model, the spatial 
waveform of the SAR target reflection can be expressed as 

G (a) = A 1 / 2 (a) X (a) (17) 

where a denotes the spatial coordinates of the along-track dimension and A (a) 
denotes the backscatter coefficient, which is also known as the reflectivity 
density. A (a) is assumed to be slowly varying. X(a) is a complex random 
process resulting from the phase variation of the numerous small sc atterers; 
it has dimensions comparable to the wavelength of the SAR. X{a) is assumed to 
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be a wide-sense stationary zero-mean complex Gaussian random process with 
independent quadrature components. The second and fourth moments of X are 

e[X(a 1 ) X* (a 2 ) ] = 8^ - a 2 ) (18) 

e[X(aj_) X* (a 2 ) X(a 3 ) X*(a 4 )] = 5(3! - a 2 ) 5(a 3 - a 4 ) 

+ 5(a! - a 4 ) 8(a 2 - a 3 ) (19) 

where 5(.) denotes a finite-impulse response which can be approximated by the 
Dirac’s Delta function when the SAR resolution width is much greater than the 
wavelength . 

C. PDF OF THE ENERGY OF A SINGLE-LOOK IMAGE 

In general, a stable sensor flight path is required during SAR mapping 
to maintain a stable sensor-target geometry and a slow drift of the Doppler 
centroid. This SAR requirement is meant to make the image processing a much 
easier task. In this case, the SAR image can be expressed as the convolution 
of the target waveform and the impulse response. The image generated from the 
i-th look, with error Af in the processing Doppler centroid, is given by 

G-(a,Af) = G (a) @ h^a^f) (20) 

It is desirable to express the above relationship as a function of time by 
assuming that the sensor is moving parallel to the target surface and with a 
constant velocity, that is 


3-5 



Gi<t,Af) = G(t) @ h ± (t,Af) 


( 21 ) 


The total energy contained in this image is given by 


r 


E 


(t, Af ) Gi*(t,Af) dt 


( 22 ) 


Notice that for simplicity, the integration along the other dimension (range) 
is not shown. For this band-limited system, E can be viewed as the summation 
of a countable number of random variables . Therefore, the probability density 
function (pdf) of Ej_ approaches a Gaussian distribution when the number of 

random variables (or the number of resolution elements contained in the image) 
becomes large. The mean and variance can be shown to be 



/• 



e [Ei (Af ) ] = 


h ± (t,Af) 

2 dt • 






/* 



v[E i (Af)] = 


hi (t, Af ) 

© hi* (t. 


Aft) dt 


A (t) dt 


( 23 ) 


( 24 ) 


The proof of the above equations is given in Appendix A. The mean of E i is 
simply the product of energy contained in the i-th impulse response and the 
integral of the reflectivity density. 

Statistically, E i and Ej are independent for i * j if the following 
conditions are satisfied for any positive integer n: 
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e [Ei (Af ) EjfAf)"] =e[E i n (Af)] etE^Af)] 


(25) 


D. PDF OF THE NORMALIZED ENERGY DIFFERENCE 


Since the arithmetic summation of Gaussian random variables is also a 
Gaussian random variable, the total energy E^ and energy difference DE also 
follow Gaussian distribution, with their mean and variance given by the 
following equations: 


e [E^ (Af ) ] 



h i (t,Af) 2 dt 


r 


A ( t ) dt 


(26) 


v [E ( Af ) ] 

T 



h L (t'A£) 


h i *(t,Af) 2 dt 


A 2 (t) dt 


(27) 



r 

[VM 

N 

„ ■ * i i ^ 

r 

e [DE (Af ) ] = 

< 

> hi(t,Af) 

| h i (tf Af) | 

► dt • 

j 


U-l 

i=N/2+l 

J 


A (t ) dt (28) 


v [DE (Af ) ] 


r 

_ hj^ (t , Af ) 

i 1 



h i *(t,Af) 2 dt 


A 2 (t ) dt 


(29) 


When the processing Doppler centroid approaches to the true Doppler 
centroid, e[DE(Af)] approaches zero. For images containing a large number of 
resolution elements, it can be proved that v [E^ (Af ) ] /e 2 [E t ( Af ) ] is also quite 
small. Under these two conditions and using the properties of the QHS given 
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above, the pdf of Ae (A f ) can be accurately modeled by a Gaussian distribution, 
with its mean and variance given by 


e [AE (Af ) ] - e[DE(Af)]/e[E (Af)] 


(30) 


v [AE (Af ) ] - v[DE(Af) ]/e 2 [£ (Af ) ] 


(31) 


Using Eqs. (15) and (16), the following results can be obtained: 


r Af 


-PBW 


„ PBW 

r — + Af 


W a (f) df - 


+ Af 


w*(f) df 


Af 


e [AE (Af ) ] 


r ^ + Af 
2 


W a (f) df 


-PBW 


+ Af 


(32) 


r 

+ Af 1 

2 




< 

w a 2 (f) df 

► • -i 


A 2 (t) dt ► 


“T + Af 





v [AE (Af ) ] = 


f 

r— + Af 1 


' 

C ' 


2 

2 




W a (f) df 

► w 4 


A ( t ) dt ► 


-PBW , A 





+ Af J 





(33) 


For small Af, the energy difference between two sides of the antenna 
pattern is nearly proportional to Af; therefore, e[AE(Af)] can be approximated 
by 
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e [AE (Af ) ] - C]_ • Af 


(34) 


where 


Ci 


2 (w a (0) - W a (PBW/2) ) 



W a (f) df 


(35) 


C x is a measure of the directivity of the antenna pattern. In Eq. (33) , 


r 

r " 




r ' 

< 

w|(f) df 

/ 

/< 


w a (f) df ► 





< J 

j 


is a measure of the time-resolution constant [13] . Hence, it can be 
represented by the product of a constant, C 2 , and 1/PBW in cases where Af is 
much smaller than PBW. Also, we may define the contrast of the target- 
reflectivity profile to be 


P = 


1 



A 2 ( t ) dt 



A ( t ) dt f 


(36) 


Thus, v[AE(Af)] can be expressed as 


v [ AE (Af ) ] = C 2 (p + 1) / (PBW • T a ) - C 2 (p + 1)/N S 


(37) 


3-9 



where T a is the time corresponding to the size of the image along the azimuth 
and N s is defined as the total number of resolution elements of the image. 


E. STANDARD DEVIATION OF THE DOPPLER CENTROID ESTIMATE 


From the statistical properties of AE(Af) given above, we can form the 
joint pdf of Af and Ae as follows: 


p(Af ,AE) 


1/DF 

C 0 ( 2jiC 2 (p + 1)/Ng) 1/2 


exp 


f - (AE - C^Af) 2 ] 
\ 2C 2 (p + 1)/N S J 


for | Af | < DF/2, | AE | < 1 


(38) 


where DF is chosen such that Eq. (37) is satisfied. C 0 is determined from the 
condition 11 p(Af,AE) • dAf • dAE = 1. The contour plot of this pdf is shown 
in Figure 3. The Doppler centroid is determined with the criterion AE = 0 . 
Therefore, the distribution of error of the Doppler centroid estimate is given 
by the conditional pdf p(Af/AE = 0) where 


p (Af /AE =0) = p (Af , A e) / p (A e = 0) (39) 

An analytical expression for p (AE = 0) is generally difficult to obtain. 
However, when N s is large enough such that p(DF/2, AE = 0) approaches zero, we 

may find that C 0 = 1 and 

p(AE = 0) = 1/ (DF • c x ) < 40 > 
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Figure 3- Contours and Projections of p(Af,AE) 



F. AZIMUTH AMBIGUITY EFFECT 

Another possible source of error in the Doppler centroid estimate is due 
to the azimuth ambiguity, which is caused by aliasing resulting from the 
limited pulse repetition frequency (PRF) employed by most radar [14] . The 
effect of the azimuth ambiguity on an SAR image is that a "ghost image" of the 
terrain outside of the beamwidth of the antenna mainlobe would be superimposed 
onto the image of the terrain within the mainlobe. 

By including the azimuth ambiguity effect, the mean and variance of the 
normalized energy difference are given by 

f ff 3 00 C 

j w a (f) df • A (t) dt - w a (f) df • A(t) dt 

k ^fi(k) Jf 2 (k) ^T rk 

e [AE (Af ) ] “ * 

'f f 3 (k> f 

< W a (f) df • A (t) dt 

k W fl ( k ) J Tk 

( 42 ) 
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v [AE ( Af ) ] 


s 


r 

-f 3 (k) 

/■ 

< 

w 2 (f) df • 

A 2 (t) dt ► 


f i (k) J 

T k J 


E 


r 


1 

< 

W a (f) df • 

A (t) dt " 


f,(k) J 

T k J 


(43) 


where 


f 1 (k) = Af + kPRF - PBW/2 
f 2 (k) = Af + kPRF 

f 3 (k) = Af + kPRF + PBW/2 (44) 

The terms T k , T£ k , and T rk are the domains of integration for targets which 
are associated with Doppler frequencies between f x (k) and f 3 (k), f x (k) and 
f 2 <k) , and f 2 (k) and f 3 (k) , respectively. 

The mean of the normalized energy difference is not very sensitive to 
the error of the processing Doppler because of the ambiguity effect. There is 
very little change to the magnitude of v[AE(Af)]. This effect leads to an 
increase in the standard deviation of the Doppler centroid estimate. However, 
the azimuth ambiguity is generally reduced by selecting a reasonably large 
PRF. In such cases, its effect on the Doppler centroid estimate is not 
significant and can be neglected. 
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G. 


SYSTEM-NOISE EFFECT ON DOPPLER ESTIMATION 


The previous sections did not treat the effect of SAR system noise on 
the DCE process. There are several types of noise sources in the SAR system: 
noise from the atmosphere, noise contributed by the receiver, and the 
quantization noise. If we assume that these noise sources are additive, the 
single-look SAR image can be expressed as 


G-(t) - G (t) @ h i (t) + n (t) 


(45) 


The bandwidth of n(t) is limited by the bandwidth of the correlator. Assume 
that n(t) is a Gaussian complex process with independent quadrature components 
and is uncorrelated with X(t) . The mean and variance of the normalized energy 
difference can be found to be 


e [AE (Af ) ] = C^Af/d + 1/SNR) 


(46) 


v [ AE (Af ) ] 


C 2 (p + 1) + 2/SNR + 1/SNR 2 
N s (l + 1/SNR) 2 


(47) 


The standard deviation of the Doppler centroid estimate is given by 
[ (C 2 (p + 1) + 2 SNR -1 + SNR -2 ) /N s ] 1/2 /C 1 . The error of the Doppler centroid 
estimate decreases when any of the following occur: 

(1) the number of resolution elements (N s ) increases 

(2) the signal-to-noi se ratio (SNR) increases 

(3) the antenna directivity (C x ) increases 

(4) the reflectivity contrast (p) decreases 
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SECTION IV 


EXPERIMENTAL RESULTS: GENERALIZED DCE ALGORITHM 


The generalized DCE algorithm, implemented with four looks, was applied 
to several Seasat and SIR-B data sets. The standard deviation of the 
estimated Doppler centroid was calculated by using the same technique 
originally described by Li [1] . A SAR data set with a size capable of 
producing a lk (range) x 4k (azimuth) pixel (full-resolution, slightly 
oversampled) image is processed to obtain 16 Doppler centroid estimates which 
correspond to spatially consecutive image blocks having 64k (range) x 4k 
(azimuth) pixels. These Doppler estimates are fitted by a linear function of 
the range. The root-mean-square deviation from the linear fit is taken to be 
the standard deviation of the Doppler centroid estimate. This method is valid 
because the actual Doppler centroid can be accurately approximated by a linear 
function of the range for both the Seasat and the SIR-B cases. The standard 
deviations of the Doppler centroid estimates from several data sets are 
tabulated in Table I. The predicted standard deviations of the Doppler 
centroid estimates and the contrasts of the reflectivity profile of the source 
are also listed in Table I. 

These reflectivity-contrast values are measured from the full-resolution 
image data by using the method described below. It is well known that the 
observed intensity of each full-resolution image pixel is the product of the 
local reflectivity of the source and a random speckle noise with an exponen- 
tial distribution. If we assume that the reflectivity of the source and the 
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Table I. Results of DCE Algorithm III and Prediction 


DATA SET 

SCENE CONTENT 

Gf dc 

p + 1 

ESTIMATION 

PREDICTION 

SEAS AT REV. 1183 
GULF OF CALIFORNIA 

UNIFORM ROUGH 
OCEAN SURFACE 

1.5 (Hz) 

1.4 

1.02 

SEASAT REV. 762 

PERIODIC 
OCEAN WAVES 

2.0 

1.8 

1.6 

SEAS AT REV. 1254 
GARDEN CITY, KANSAS 

AGRICULTURE 

FIELD 

2.7 

2.4 

3.0 

SEASAT REV. 351 
SAN GABRIEL MTNS . 

MOUNTAINS 

4.1 

2.8 

4.2 

SEASAT REV. 351 
L.A. AIRPORT, CALIF. 

AIRPORT IN 
URBAN AREA 

8.1 

3.7 

■ 

SIR-B REV. 6 
W. MONTREAL, QUEBEC 

URBAN AREA 

10 f 4 

4.6 

11.2 
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random speckle noise are independent, the contrast of the source reflectivity 
is given by 



where M x is the mean pixel intensity and M 2 is the mean squared pixel 
intensity. The antenna directivity constant and the effective bandwidth 
constant C 2 are calculated based on the given antenna pattern [15] . 

Table I shows that the predicted Doppler centroids are very close to the 
Doppler centroids estimated from the SAR data for a uniform ocean surface, 
periodic ocean waves, and agriculture fields. This indicates that the 
performance analysis given in Section II does apply well to the QHS . Table I 
also shows that large differences between predictions and estimates exist for 
the mountain and urban sources. This indicates that these sources are not 
accurately modeled by a QHS. However, the prediction does provide an 
indication of the growth of estimation errors for these cases. 
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SECTION V 


OPTIMAL DOPPLER CENTROID ESTIMATION FOR QHS 


The DCE algorithms described above do not yield an optimal Doppler 
centroid estimate for a QHS. This fact is illustrated by the following 
example. Consider a source consisting of two equal areas, each of which has a 
uniform reflectivity and a very large signal-to-noise ratio. Assume that the 
reflectivity level of one area is much greater than the other. The errors of 
the Doppler centroid estimates obtained from each individual area are almost 
equal; therefore, arithmetic averaging reduces the estimation error by a 
factor of I/V2. However, the error of the Doppler centroid estimate obtained 
by applying those algorithms described in Section IIB over the entire area 
would be almost equal to those from each individual area. 


In the above example, the greater uncertainty contained in the Doppler 
centroid estimate from those algorithms is due to the fact that the weight is 
given by the percentage of energy contained in each subarea instead of by the 
size of each subarea. This can be seen from the following expression for AE : 


E T (1) _ E t (2) 

AE = AE(1) + AE (2 ) 

tl rn Ej rn 


(49) 


where AE(1) and AE(2) stand for the normalized energies corresponding to those 
subareas, and E T (1) and E x (2) are the total energies corresponding to those 

subareas . 
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The above example also gives a hint for an optimal approach. A QHS can 

be thought of as a source consisting of a set of spatially segregated 

subsources, each of which is a homogeneous source. The first step of an 

optimal DCE process for a QHS {also called an optimal QHS-DCE process) is to 

obtain a set of optimal Doppler centroid estimates from these subsources. The 

next step is to obtain a final estimate based on these preliminary results. 

The solution for the second step is quite obvious. An optimal Doppler 

centroid estimate based on the maximum likelihood criterion is given by the 

weighted sum of those preliminary estimates, where each weight coefficient is 

inversely proportional to the variance of each preliminary estimate. Let 

a- <k) be the standard deviation of the Doppler centroid estimate of the k-th 
r dc 

subsource. The associated weight coefficient to be used in the optimal QHS- 
DCE process is given by 


W c (k) 



00 


I 



(k) 


( 50 ) 


The standard deviation of the final Doppler centroid estimate is given by 


f K 


<*f 


dc 


YL/2 


W 2 (k) CTf (k) 

^ '‘■dc 


k=l 


f K 


yi/2 


CTf 2 (k) 

■ L dc 


k=l 


( 51 ) 
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The optimal homogeneous-source DCE algorithm using a linear estimation 
approach has been developed and will be reported in another publication [16] . 
The detailed analysis of this algorithm will not be repeated here. However, 
an introduction to the problem and the optimization analysis involved are 
given in Appendix B. In the following we describe the processing procedures, 
which include both time domain and frequency domain approaches. The time 
domain approach can be incorporated in the optimal QHS-DCE process more easily 
than the frequency domain approach. 

A. OPTIMAL DOPPLER CENTROID ESTIMATION FOR HS 


The frequency domain approach is used to obtain the range-averaged power 
density spectrum from a full-resolution image. The next step is to convolve 
this spectrum with W a (f)/W a (f), where W^(f) is the derivative of W a (f) with 
respect to f. After convolution, the frequency associated with a sample, 
having its value closest to zero, is the Doppler centroid estimate. 

The time domain approach is to process the SAR echo response into N 
independent single-look images, each with 1/N of the total processing 
bandwidth, where N is an even number. These single-look images are weighted 
in the frequency domain with W a (f)/W a (f), by incorporating this weight 
function in the reference function. Compute AE according to Eq. (4) . The 
Doppler centroid estimate is then obtained by subtracting the estimated 
processing Doppler error Af (given in Eq. (52)) from the processing Doppler. 
Repeat this procedure until AE approaches zero. The estimated processing 
Doppler error is given by 
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PBW 

2 


w;<f> 


-PBW W a ^ 


df 


Afdc = 


PBW 

2 fw; ( f)Y 


-PBW 


V W a(f) j 


df 


AE 


(52) 


The standard deviation of the Doppler centroid estimate obtained from 
this algorithm is given by 


PBW 

2 (K^i 

1 df /PBW 


Vl/2 


-PBW 


v w a ( f ) 


a f d c " 


f (SNR) 


(53) 


N 


1/2 


The polynomial function f (SNR) is determined by the numerical analysis method. 
For a large SNR, the value of f(SNR) is equal to one. 


In the implementation of the optimal QHS-DCE algorithm, determination of 
subsource boundaries is not necessary. A simple method is to use the 
intensity histogram of the multi-look image. The procedure includes the 
following: Perform multi-look summation to generate the N-look image (to 

reduce speckle uncertainty) and its histogram. Divide the intensity range, 
covering more than 99% of the total image pixels, into K intensity intervals. 
Preliminary Doppler centroid estimates are then obtained from subsets of 
pixels having the same intensity level, A detailed block diagram of this 
optimal QHS-DCE process is illustrated in Figure 4. 
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Figure 4. Block Diagram of an Optimal QHS-DCE Algorithm 











SECTION VI 


EXPERIMENTAL RESULTS: OPTIMAL QHS-DCE ALGORITHM 

The optimal QHS-DCE algorithm was applied to the same SAR data sets 
listed in Table I. The same linear-fit method was used to estimate the 
standard deviation of the Doppler centroid estimate. Results are listed in 
Table II. In these cases, the center portions of the linear-fitted lines 
obtained from both algorithms differ by an amount much smaller than the 
standard deviation listed in Table II. An improvement of the estimation 
accuracy can be seen in all the cases, even in urban and mountain sources. 

However, for urban and mountain sources, the Doppler centroid estimates 
are not as accurate as those obtained from the ocean and agriculture areas. 
This could be due to the fact that the reflectivity profiles of these sources 
contain high-frequency variations, which are beyond the assumptions of a QHS . 
Also, the assumption of phase independency among numerous small scatterers may 
not be true for the urban sources. It is believed that these factors degrade 
the accuracy of the Doppler centroid estimates. 
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Table II. Results of an Optimal QHS-DCE Algorithm 


DATA SET 

SCENE CONTENT 

Gf dc 


ALGO. 

III 

OPTIMAL 

SEAS AT REV. 1183 
GULF OF CALIFORNIA 

UNIFORM ROUGH 
OCEAN SURFACE 

1.5 

(Hz) 

1.0 

SEASAT REV. 762 

PERIODIC 
OCEAN WAVES 

2.0 

1.2 

SEASAT REV. 1254 
GARDEN CITY, KANSAS 

AGRICULTURE 

FIELD 

2.7 

1.5 

SEASAT REV. 351 

SAN GABRIEL MOUNTAINS 

MOUNTAINS 

4 . 1 

2.8 

SEASAT REV. 351 
L. A. AIRPORT, CALIF. 

AIRPORT IN 
URBAN AREA 

8.1 

3.3 

SIR-B REV. 6 
W. MONTREAL, QUEBEC 

URBAN AREA 

10.4 

5.5 
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SECTION VII 


CONCLUSION 


We have presented a DCE algorithm which is generalized from previously 
reported algorithms. Detailed performance analysis was given and validated by 
the experimental results. This performance analysis is directly applicable to 
Li f s algorithm. In Curlander's algorithm, the Doppler centroid estimate can 
be slightly more accurate when the selected processing Doppler values cover a 
wide range (to increase the effective total number of resolution elements) and 
the azimuth ambiguity levels contributed from two sides of the main targets 
still balance well. 

In the second part of this publication, an optimal DCE algorithm for the 
quasi-homogeneous source was presented. Experimental results show that it 
does improve the accuracy of the Doppler centroid estimates for all the tested 
data sets . 
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APPENDIX A 


PROOF OF (23) AND PROOF OF (24) 


A. PROOF OF (23) 


From Eqs . (21) and (22), the mean of E^Af) can be expressed by 


r r 


e (Ei (Af ) ] = 


e[G(t')G* (t") ]hi(t - t',Af) 


h?(t - t" , Af ) dtdt'dt" 


(Al) 


By use of Eqs. (17) and (18), it can be shown that 


e [G(t 1 ) G* (t") ] = A 1/2 (t')A 1/2 (t")8(t’ - t") 


(A2) 


Substituting (A2) in (Al) , we have 


e [Ej^ (Af ) ] = 


h^t - t ' , Af ) A(t')dt'dt 


(A3) 


In (A3), h^ (t - t',Af) can be viewed as a delta function because of the 
assumption that A(t') is a slowly varying function. This implies 


e[E i (Af)] 


hj^tfAf) 2 dt • 


A (t) dt 


( A4 ) 
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B. 


PROOF OF (24) 


From Eqs . (21) and (22), the mean of E? (Af ) can be expressed by 


e[Ef (Af) ] 


r 


r 


x 


• • • • e[G(t 1 )G*(t 2 )G(t 3 )G*(t 4 ) ] 


h L (t - t x , Af ) h| (t - t 2 , Af ) 


h^t' - t 3 , Af ) tq (t ' - t 4 , Af ) 


dtdt 1 dt 1 dt 2 dt3dt 4 


(A5) 


By use of (17) and (18), it can be shown that 

1/2 1/2 1/2 1/2 

e[G(t 1 )G* (t 2 )G(t 3 )G*(t 4 ) ] = A (t x )A (t 2 ) A (t 3 )A (t 4 ) 


[6 { 1 3 - t 2 )5(t 3 - t 4 ) + 5(t 3 - t 4 )5(t 2 - t 3 ) ] (A6) 


Substituting (A6) in (A5) , we have 
e [E? (Af ) ] = e 2 [E i (Af)] 


r r 


+ 


*i<ti 


j j 


t 2 ,Af) @ h* (t x - t 2 ,Af) 


A (t x ) A (t 2 ) dt x dt 2 

(A7) 


A-2 


1 1 ■mil Mil 


In (A7 ) , h i (x) h? can be viewed as a delta function because of the assumption 


that A (t f ) is a slowly varying function. This implies 


v[Ei(Af) ] 


/• 

i i 2 

I hi (t, Af ) @ i (t, Af ) I dt 


A 2 (t)dt 


J 


(A8) 
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APPENDIX B 


OPTIMAL DOPPLER CENTROID ESTIMATION FOR HS 

The spectrum of the SAR response from a homogeneous source can be 
expressed by the following function: 

S (f ) = pW(f - f dc ) X(f)X*(f) (Bl) 

A 

where p is proportional to the power of the response and X(f) is the Fourier 

a 

transform of random process X(t) in Eq. (17) . X(f) can be shown to be a zero- 

mean complex Gaussian random process with independent quadrature components . 

/\ 

The second and fourth moments of X are of the same form as those of Eqs . (18) 

A A 

and (19). W(f - f dc ) is equal to the antenna pattern W a (f) given in Eq. (9). 

The Doppler centroid estimate from a linear estimator can be expressed 
by 


f dc< f l> = f l + Af < f l> (B2> 


where 


Af (f x ) 


r f 1 +F/2 

S (f ) g (f - f x ) df 

Jfj-F/2 


(B3) 
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Here, g(f) is the weight function; f x is the frequency reference of g(f). F 
is the frequency bandwidth over which the estimation is performed. The mean 
and variance of Af (f^ can be expressed by the following equations when f x 
approaches f dc . 


e [Af (f x ) ] » P 


F/2 


W(f)g(f) df + <fj. - f dc ) • P 


-F/2 


F/2 


-F/2 


W (f)g(f) df (B4) 


and 


v[Af (f x ) ] 


C F/2 


W 2 (f ) g 2 (f ) df 


-F/2 


(B5) 


Here, we shall find a function for g(f) which yields an unbiased and 
optimal estimate of f dc . The "optimal" estimate is interpreted here as the 
one having the least mean-square error. Therefore, g(f) must minimize 
vtAf^)] and satisfy the following: 

" F/2 

W(f)g(f) df = 0 1 

J -F/2 


'F/2 


w' (f)g(f) df = 1 


-F/2 


(B7) 
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Let g(f) = C g x (f) g 2 (f)r and solve for C from Eq. (B7) . Substituting C in 


Eq. (B5) , the problem becomes one of maximizing 
fF/2 

" W' (f) gi (f)g 2 (f) df" 

* -F/2 


" F/2 

W 2 (f)g2(f)g2(f) d f 

J -F/2 


(B8) 


By choosing g^(f) to be equal to 1/W(f) and using the Cauchy-Schwarz 
Inequality, we can show that the solution for g 2 (f) is the function 

A f A 

W (f)/W(f) . Hence the weight function g(f) is given by 


g(f) 


1 


F/2 


-F/2 



w' (f) 
w 2 (f) 


(B9) 


A 

For a conventional SAR antenna pattern, W(.) is symmetric and its first- 
order derivative is asymmetric, i.e.. 


w (f) - - w (-f) 


(BIO) 


This implies that g(f) is also asymmetric and satisfies (B6) . 
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